Als Typoskript gedruckt 

Technische Universitat Dresden 

Herausgeber: Der Rektor 



(N 

o 

(N 
u 

< 

(N 



CM 

< 



The Elusive Drude-Born-Fedorov 

> 

o ; Model for Chiral Electromagnetic 

£j : Media. 

q \ Henrik Freymond and Rainer Picard 

; Institut fur Analysis 

MATH-AN-03-2012 

X 



The Elusive Drude-Born-Fedorov 
Model for Chiral Electromagnetic 

Media. 



Henrik Freymond & Rainer Picard 
Institut fur Analysis, Fachrichtung Mathematik 
Technische Universitat Dresden 
Germany 
rainer.picard@tu-dresden.de 



Electro-magnetic wave propagation in more complex linear materials such 
as bi-anisotropic media have come to a considerable attention within the last 
fifteen to twenty years. The Drude-Born-Fedorov model has been extensively 
studied mostly in the time-harmonic case as a model for chiral media. In the 
physically relevant time- dependent case the record is much less convincing. In 
this paper we focus on this case and analyze the Drude-Born-Fedorov model in 
the light of recently developed Hilbert space approach to evolutionary problems. 
The solution theory will be developed in the framework of extrapolation spaces 
(Sobolev lattices). 

Keywords and phrases: Maxwell's equations, bi-anisotropic material, chiral media, extrapolation 
spaces, Sobolev lattices 
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Introduction 

The evolutionary system of the macroscopic Maxwell's equations 

curl H — 3qD = j , 
curl E + d B = 0, 

complemented by suitable boundary conditions such as the vanishing of the tangential 
component of the electric field E and initial data for D and B is completed by a material 
law connecting the combined six component vector field (D, B) composed from the electric 
displacement current D and the magnetic induction B with the electro-magnetic field 
(E,H) (do denotes the time derivative). In the simplest case it is assumed that 

D = eE 
B = u,H 

with constant coefficient £,|ie M >0 (permittivity e, permeability /i). This case is referred 
to as the case of isotropic and homogeneous media. There is a long history of various 
description of more complex electromagnetic media and a shorter one for their mathemat- 
ical treatment. As - in a sense - the conclusion of the treatment of the anisotropic and 
inhomogeneous media case we refer to the presentation of the Maxwell system in the book 
by R. Leis [9j from 1986 and the references given there. The functional analytical approach 
utilized shows that indeed a mathematical solution theory can be obtained for rather gen- 
eral media merely requiring that e,fi are time-independent, bounded, self-adjoint, strictly 
positive mappings in L 2 -type spaces of vector fields in an non-empty open set Q c ]R 3 
containing the media characterized by e and /i, which covers the case of real multiplicative 
matrix-valued operations with in Q L°°-bounded entries as a special case. Although, the 
necessary concepts are - as one says - well-known since the late seventies of last century, it 
seems not too widely adopted that by appropriately generalizing the concept of assumption 
of boundary conditions boundary singularities such as corners, edges, cusps, even fractal 
boundaries can be included. The price to be paid for this is a more subtle approach ex- 
ploiting the full power of functional analytic concepts (as a general reference for functional 
analytical concepts see [23]). As it turns out, however, this can be achieved purely 
in a Hilbert space setting, thus reducing conveniently the conceptual complexity of our 
considerations. 

Writing ( • | • ) for an L 2 -type inner product regardless of the number of components, we 
have 



E 
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for vector fields 




E = (Ei, E 2 , E 3 ) = E 2 \, H = (Hi, H 2 , H 




m L 2 (S}). For ease of notation we shall rarely note the number of components and write 
e.g. L 2 (St) instead of L 2 (SI) 3 , since the number of components will always be clear from 
the context. It is clear from the arguments e.g. in [13], although the impact and relevance 
of this observation for the more recent discussion of electromagnetic meta-materials could 
not be realized at the time, that the mathematical methods developed in connection with 
anisotropic, inhomogeneous media extend to the case of material laws of the block matrix 
form 

with 

in an obvious way, provided that the entries e, [i, (5 are time-independent, bounded, linear 
mappings in L 2 -type spaces of vector fields such that the block matrix operator 

M : L 2 (SI) -»■ L 2 (SI) 

E\ ( e P\ ( E 



H J \P* ft J \H 
is still self-adjoint and strictly positive!]] 

Based on a structural observation described in [T7] we shall here consider in the following 
a class of well-posed evolutionary equations in even more complex media. Such media 
have become of strong interest in recent decades due to their prospects in view of their 
interesting physical properties (for a survey we refer the reader to [TT| [8]). This problem 
class, which conveniently contains materials with memory, will lead us to the so-called 
Drude-Born-Fedorov model for chiral media, which has been initially introduced as an 
approximation of the actual material description in an attempt to avoid non-local-in-time 
terms (optical response). 

The model starts from a material relation of the form 

D = e(E + r] curl E) , B = fx(H + r) curl H) , 
where e, /i 6 M>o and rj e R is the so-called chirality parameter. 

There is a long list of publications dealing with this model. Despite this fact and contrary to 
various claims made, it appears that this model has not been fully understood in sufficient 



1 This is the case if and only if e, fi — (3*e 1 (3 are selfadjoint and strictly positive, which amounts to a 
smallness condition on the off-diagonal entry f3. 
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1 Space-Time Evolution Equations 

generality. It is due to this state of affairs that we have chosen to speak of the "elusive" 
Drude-Born-Fedorov model. So has it only more recently been understood that the model 
has to be considered in the - much less forgiving - dynamic case, O [221 E El E], rather 
than the so-called time-harmonic case. 

We will show that in the framework of our approach from [T7] the Drude-Born-Fedorov 
model will have to be considered as a rather degenerate case, which, however, by some 
additional considerations can be elegantly embedded into the standard theory. Indeed, 
we are led to a simple evolution equation with a bounded generator. This has in essence 
already been observed in [10] for simply connected domains with smooth boundary. It has 
gone unnoticed, however, that the result can be extended to fairly arbitrary open sets with 
a boundary permitting a suitable local compact embedding result, [T^l IT5| [5]. As a by- 
product we can discuss well-posedness for a class of bi-anisotropic materials with modified 
Drude-Born-Fedorov type material behavior, some cases of which have found previous 
attention, compare [T2J Morro 2002], [211 Sjoberg 2008]. 



1 Space-Time Evolution Equations 

1.1 A brief summary of Sobolev chains and lattices 

We recall first the standard construction of chains of Hilbert spaces associated with a 
normal operator O in a Hilbert space X with in its resolvent set, see e.g. [IB]. Defining 
H k (O) as the completion of D (O fc ) with respect to the norm | • \ k induced by the inner 
product ( • | • )fc gi ven by 

(f,g) ~> (O k f\O k g) x 

for k e Z, we obtain a chain (H k (0)) keZ of Hilbert space, where the order is given by 
canonical, continuous and dense embeddings 

H k+1 (O) — H k (0),ke Z. 

Note that Hi (O) is the domain of O and for k e N 

H k (O)^X^ H_ k (O) 

is a Gelfand triple. By construction 

H x (O) - H (O) 
f~Of 

is a unitary mapping and it follows by induction that for k e N also 

H k+1 (O) - H k (O) 
f~Of 
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1.1 A brief summary of Sobolev chains and lattices 



is unitary. Moreover, O extends to a unitary mapping 

H k+1 (O) - H k (O) 
f~Of 

also for negative k and so to all fceZ. 

Remark 1.1. Since = d\zO + \ 3m0 is a normal operator in H (O) we have 

O* = ^HeO-iamO 

Since D (O) = D (O*) and 

|0/|* = \O m f\ x 

that also O has a unitary extension 

tf fe+a (O) - fr fc (O) 

for k e Z. Keeping the notation O* for this extension, as we shall do, leads to a possible 
confusion, since O* could be misinterpreted as the adjoint, i.e. the inverse of O as a 
unitary mapping. To avoid this misunderstanding we shall not write the inverse of a 
unitary operator as its adjoint. 

With this convention we can safely say that we have established s and (0*) s as unitary 
mappings 

H k+S (O) - H k (O) 
H k+S (O) - H k (O) 

for s,k s Z. 

We shall utilize this abstract construction for various special cases. Note that for O e 
L (H, H) we have 

H s (O) = H 

as topological linear spaces with merely different inner products (inducing equivalent 
norms) in the different Hilbert spaces H s (0), s e Z. This indicates that considering 
continuous linear operators O does not lead to interesting chains. 

A particular instance of this construction is if O is the time-derivative do. We recall, e.g. 
from |17[ 118] . that differentiation considered in the complex Hilbert space H g o(M) := {/ e 
L^ oc (R)\(x h> exp(— gx)f(x)) e L 2 (R)}, g e M\{0}, with inner product 

(f,g) (f,g) e ,o : = f(x)*g(x) exp(~2gx) dx 



and 
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1 Space-Time Evolution Equations 



can indeed be established as a normal operator, which we denote by do i6 , with 

^ do, e = Q- 

For 3m do iB we have as a spectral representation the Fourier-Laplace transform £ g : H Sj o(M) - 
L 2 (R) given by the unitary extension of 

Coo(M) <=tf ei0 (K) ^L 2 (R) 

^ •-»■ J exp (-ixy) exp (-£>y) 0(y) dy^j . 

In other words, we have the unitary equivalence 

3md 0i g = C~ l mC e , 

where m denotes the selfadjoint multiplication-by-argument operator in L 2 (M). Since 
is in the resolvent set of 3q we have that 8q is an element of the Banach space 
L (H ey o(M.), H g fi(M>)) of continuous (left-total) linear mappings in H Bj o(M.). Denoting gen- 
erally the operator norm of the Banach space L (X, Y) by || • \^xy)^ we § e ^ f° r ^cT 1 

ll^ll =1 
Not too surprisingly, we find for g > 

fa»(*) = f <p(t) dt 

J— co 

and for g < 

roo 



for all 99 e (M) and jeK. Since we are interested in the forward causal situation, we 
assume g > throughout. Moreover, in the following we shall mostly write 8q for 8q^ if 
the choice of g is clear from the context. 

Thus, with O = do, e in the above general construction we obtain a chain (H e ^ (K)) fcgZ of 
Hilbert spaces 

H e> k QR) := H k (d 0>e ) . 

Similarly, with O = im + g in L 2 (K) we obtain the chain of polynomially weighted L 2 
spaces 

with 

L\ (R) := {/ e L 2 ' loc (E) I (im + / e L 2 (R)} = ff fe (im + e ) 

for fceZ. 
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1.1 A brief summary of Sobolev chains and lattices 



Since the unitarily equivalent operators <9o, e an d im + q can canonically be lifted to the 
X-valued case, X an arbitrary complex Hilbert space, we are lead to a corresponding chain 
(Hg^k (R, X)) keZ and [L\ (R, X)) fegZ of X-valued generalized functions. The Fourier-Laplace 
transform can also be lifted to the X-valued case yielding 



as a unitary mapping for fceN and by continuous extension, keeping the notation C Q for 
the extension, also for k e Z. Since £ e has been constructed from a spectral representation 
of 3mdo tB , we can utilize the corresponding operator function calculus for functions of 
3mdo t g. Noting that do tB = i3mdo tB + g is a function of 3mdo tB we can define operator- 
valued functions of do- 

Definition 1.2. Let r > > and M : Bc(r,r) — > L(H,H) be bounded and analytic, 
if a Hilbert space. Then define 



for 06 (R, 

Remark 1.3. The definition of M^q -1 ) is largely independent of the choice of g in the sense 
that the operators for two different parameters Q\ , g 2 coincide on the intersection of the 
respective domains. 

Simple examples are polynomials in d^ 1 with operator coefficients. A more exotic example 
of an analytic and bounded function of d^ 1 is the delay operator, which itself is a special 
case of the time translation: 

Example 1.4. Let r > 0, g > hsR and u e H es o(M., X). We define 

Th.u := it( • + h). 

The operator e L(H gt0 (M., X), H b q(M., X)) is called a time-translation operator. If h < 
the operator is also called a delay operator. In the latter case the function 

Bc(r,r) 3 z >— > M(z) := exp(,z _1 /i) 

is analytic and uniformly bounded for every r e R >0 (considered as an L (X, X)-valued 
function). An easy computation shows for u e H B) q (R, X) that 

u( ■ + h) = L*exp((im + g)h)~h B u = M(5 ( 7 1 )w = exp((5 _1 ) 1 h) u. 
This shows that 



H e>k {R,X) - (R,X) 
f~£ e f 




where 




r ft = exp (h/d x ) = exp (ftd ) ■ 
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1 Space-Time Evolution Equations 



Another class of interesting bounded analytic functions of d are mappings produced by 
a temporal convolution with a suitable operator-valued integral kernel. 

Let now O denote a normal operator in Hilbert space H with in the resolvent set. Then 
O has a canonical extension to the time-dependent case i.e. to H e Q (R, H). Then do and O 
become commuting normal operators and by combining the two chains we obtain a Sobolev 
lattice in the sense of |19] based on (do, O) yielding a family of Hilbert spaces (g e R\ {0}) 

(H eik (R, H s (0))), eZ 



\d k O s v\ 



with norms | • \ ks given by 

v 

for k, s g Z. Note that 

H (O) = H 
independent of the particular choice of O. 



1.2 Abstract initial value problems 

We shall discuss equations of the form 

(d M (do 1 ) + A)U = J + 5®W , (1) 

where for simplicity we shall assume that A is skew-self adjoint in a Hilbert space H and 
M is a regular material law in the sense of [T6| IT9] . More specifically we assume that M 
is of the form 

M (z) = M + zM 1 (z) 

where Mi is an analytic and bounded L (H, H)-v&lued function in a ball (r, r) for some 
r e M >0 and M is a continuous, selfadjoint and strictly positive definite operator in H. 
The operator M (d^ 1 ) is then to be understood in the sense of the operator- valued function 
calculus associated with the selfadjoint operator 3m (do) = ^ (do + d^). 

For the data we assume 

J 6 H e>0 (R, H) , J = on R <0 , 

and 

W eH, 

which makes (pQ) an abstract initial value problem. The appropriate setting turns out to 
be the Sobolev lattice 

\H e , k (r, H s (^M^A^JmY 1 + l) ) ) 

where, however, only the spaces with s = — 1, 0, 1 and k = —2, —1,0, 1 are actually utilized. 
From [TB| [T5] we paraphrase the following solution result on which our approach to the 
Drude-Born-Fedorov model can conveniently be based. 
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1.2 Abstract initial value problems 



Theorem 1.5. The abstract initial value problem (T7J) has a unique solution U e H g _i (R, H). 

Moreover, 

F~(d M (d^ 1 )+A)- 1 F 

is a linear mapping in L (H Q ^ (R, H) , H e k (R, H)) , k e Z. These mappings are causal in 
the sense that if F e (R, H) vanishes on the time interval ] — go, a], then so does 
(d M (do 1 ) + A) 1 F, a e R, k e Z. In particular, we have 

*7 = on R <0 . 

To link up with a more classical interpretation of the assumption of the initial data we 
record the following regularity result. 

Theorem 1.6. Let U e H e -\ (R, H) be the unique solution of the abstract initial value 
problem (QJj. Then, we also have 

UeH e>0 (R,H) (2) 
and U is a continuou^ H -valued function on R\{0}. Moreover, 

U (0+) = Mo" 1 Wo, (3) 

where the limit is taken in H. 

Proof. These stronger regularity statements will rely on the function calculus for the skew- 
selfadjoint operator sj M^Ay/ M^ 1 or - depending on the point of view - on one-parameter 
semi-group arguments. Let U s H e -\ (M.,H) be the solution of ([!]). Noting that 



(d M + Ay 1 = Jm^ 1 d + Jm^aJm, 1 Jm ~ 1 



we have 



(d + x/m^M, (do 1 ) ^JMo 1 + x/m^A^/m^ 1 ^ VMoC/ = \[m^3 + S ® a/m^Wq. 



From this we have 



^1 + (d + x/m^A^/m^ Mi (do 1 ) \ y/MoU = (do + yf^Ayf^j \Jm~ x J + 



+ {do + JM^AJMq 1 ) 5 ® \/M, 



-i 



-l 



2 In the usual sense of having a continuous representer. 
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1 Space-Time Evolution Equations 



We have 



U + J^A^mA S^x/m^Wo = (t - x«, (*) exp (-tyJ^AyJ^j y^Wc 

by uniqueness of solution, which as we read off is in H g (R, i/) and is a continuous function 
on M\{0}. In particular, we also have 



(d + y[^Ay[M^j 6 ® y^WoJ (0+) 



(4) 



-l 



in H. Moreover, for g e IR>o sufficiently large we have that ydo + y/ M 1 Ay/ M 1 J Mi (<9 x ) 
is a contraction since 

-l 



(do + y[^AyfM^\ f 



e,o,o 



e,o,o 



Therefore 



[/ = yMo- 1 (1 - Q)- 1 f d + yj Mq 1 A\J M Q l j [yJM Q l J + 5 ® yj Mq X Wq\ e H gfi (E, if) 
where Q = - (d + y/M^Ay/M^ * M x (c^ 1 ) . We also have that 

1 + (do + ^fM ^1 A x fM^^ M, {d, 1 ) \ = 

= l-(d + ^fM ^1 A x fM^^ M x {d, 1 ) (\ + (do + y[^A^M^ M, [d, 1 ]^ 



and so 



U = JMv 1 do + Jm^AJM, 1 5 (x) JM, 1 W + 



(5) 



+yjM 1 [d + yjM^AyjM, 1 ) F 
for some F e H gfi (K, H) with F = on R <0 . Since 
f exp (- (t - s) y/M^AyfM^^j F (s) ds < f \F (s)\ ds 

J— co V / J— oo 

= exp (gs) \F (s)\ exp(-gs) ds, 

J— co 



< 



|F(s)| 2 exp (— 2gs) ds, 
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existence of the integral Sf ^ exp (^— (t — s) y/ M 1 Ay/ M ^ F (s) ds and continuity of t >— > 
S/ ^ exp (^— (t — s) yj M ' 1 Ay/ Mq 1 ^ F (s) ds is clear. Since F = on lR<o we have in 



particular that 



J exp f- (t - s) x/m^Ax/m^^J F (s) ds = 0. 
Since the first term of satisfies property (jl]) we have indeed 



17(0+) = x/M^x/M^Wo = M^W . 

□ 

Remark 1.7. In the case A = even stronger regularity follows. Indeed, the solution is 
then such that 

U- Xu>0 ®M^W Q eH eA (R,H). (6) 

2 Discussion of the Drude-Born-Fedorov Model 

We begin by analyzing the Drude-Born-Fedorov material relation (e, \x 6 R>o, r] e M\ {0}) 

D = e (E + 7] curl E) , B = yu (H + 7] curl H) 

in the light of the above general theory by substituting Maxwell's equations back into the 
Drude-Born-Fedorov relation (ignoring possible source terms) to obtain a modified material 
relation of the form 

D = e (E-tj d B) , B = n (H + r\ d D) . 

This is a modified Condon model, jU Condon 1937], compared to which in the right-hand 
sides B and D are replaced by H and E, respectively. 

Although similar in concept to the Condon model, where replacing a convolution term 
by a two-term Taylor approximation leads to no reasonable results, here this type of "ap- 
proximation" can actually be justified, see [HI Frantzeskakis, Ioannidis, Roach, Stratis, 
Yannacopoulos (2003)]. With a slight reformulation we have Maxwell's equations 

and the Drude-Born-Fedorov material relation assumes the formal shape 

1 V^o\(^D\/^=E 
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2 Discussion of the Drude-Born-Fedorov Model 



or rather 



In the language of the above theory 



1 -rj v^^o 
77 V^^o 1 



1 + ( -^^V 1 



1 \ 

2\ 



Note 







*He 1 " ^ = 

\ '1 v^M / 

which shows that the resulting equation is not covered by the above theory directly and 
actually is not a differential equation in time at all: 




Moreover, this whole consideration is assuming that we consider curl to be a self-adjoint 
realization in L 2 (JV), which can only be achieved in very special cases, such as Q = M?\N, 
where N is a set of capacity zero. (e.g. in Q = R 3 ). 

However, for media occupying an arbitrary open subset Q of M 3 we fortunately have a nat- 
ural choice of boundary condition, which turns curl with a corresponding choice of domain 
into a selfadjoint operator. To formulate this condition properly we need to introduce curl 
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as the closure in L 2 (f2) of curl restricted to (f2) vector fields (we do not indicate the 
number of components). Then curl is properly defined as the adjoint of curl: 

curl := ( curl 



Similarly, div is defined as the closure in L 2 (£7) of div restricted to C x (Q) vector fields. 
Then 

grad := — (div) 

is the usual weak derivative in L 2 (fi), , see [T^ 7rk 1998] for the conceptual details. Con- 
tainment of a field E in D (curl) is the proper weak generalization of the classical boundary 

condition "n x E = Oon dQ", whereas E e D (div) generalized the classical boundary con- 
dition "n-E = Oon d{l". It is important here to keep in mind that no regularity assumptions 
on the boundary and no trace results are needed for these generalized constructions. 

A suitable boundary condition for the Drude-Born-Fedorov model can now be stated in 
terms of the range of curl, R (curl) , and of the domain of div, D (div) . 

We require 

curl E e R (curl) (7) 

or equivalently 

curl E e D (div) (8) 

and 

cmAE 1U N , (9) 
where Kjv denotes the set of harmonic Neumann fields 



H N = [Ee D (div) | div£ = 0,curl£ = oj 



Condition (jSJ) generalizes the classical boundary condition "n • curl E = on dtt" to non- 
smooth boundaries and data[§. We shall denote the operator curl subject to boundary 

o o 

conditions (JHJ), by curl. Under fairly general assumptions it can be shown that curl is 
actually selfadjoint and even has - apart from - only discrete spectrum ad I curl 



a p curl \\{0} = a d curl . (10) 



Note that assuming data in D (^divj boundary condition (|8]) is induced by the requirement that 

E e D (div^ . Moreover, in the simply connected case Hn = {0} so that ([7]) reduces to ©. These 
observations are the link to the set-up utilized in [10] . 
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2 Discussion of the Drude-Born-Fedorov Model 

The most general result, Filonov 2000], merely requires Q to be an open set with 

o 

bounded measure to obtain the same properties for curl. Exterior domains, i.e. open 

o 

sets with compact complement, also support selfadjointness of curl, [151 nk 1998], if the 
boundary satisfies a local compact embedding property, see |20j . Indeed, although in the 

o 

exterior domain case the spectrum of curl is not purely point spectrum anymore, we still 
maintain (UOp . 

In order to avoid technicalities and to immunize the results presented here against possible 
future improvements, we make the crucial property of selfadjointness our core assumption 
for the analysis of the Drude-Born-Fedorov model: 

Hypothesis f2: We assume that Q is an open subset of M 3 such that 

o 

• curl is selfadjoint. 

Under this general assumption we shall re-inspect the Drude-Born-Fedorov model from 
another perspective. We first note that 



D = e(E + r] curl E) = (1 + rj curl) eE, B = fi(H + rj curl H) = (1 + n curl) \xH. 
So, imposing our general Hypothesis Q the material relation takes on the form 

D b)-( 1 + ^ ;a )(11)(h 



The initial value problem for Maxwell's equation now reads formally 



Now let P r , denote the orthogonal projector onto R^l + n curl) . Consequently, (1 — P 7] ) 
is the orthogonal projector onto N (^1 + rj curl) and so 

(1 - P v ) curl = rj- 1 (1 - P n ) curl + 1 - l) 

= -v- 1 (1 - P v ) (11) 

curl (1 - P v ) = rj' 1 ^r/curl + 1 - l) (1 - P v ) 

= T]- 1 Ucil + l) (1 - Pr,) V' 1 (1 - Pr,) 

= -rj' 1 (1 - P„) . 



16 



Thus, P v and curl commute and we have 

I + 77 curl) d pj £ n °) \ E U ) + 



fj, J \H 



l + ,curl^oP 7? (^)p,(|) + 
+ (1 + 7, curl) ^(oj)(l-^)(f) + 



+ ( " n ) P„curl = P^J + <y® P,Wb 



° _1 )Pcurl( ^ 
° "q 1 ) (l-i^)curl( ^ ) 



and so 



— 1 \ 1 E 

+ ( 1 r, 7T„CUrl7r*7r„ rjr I = 7T„ J" + 5 (X) 7T„W 



It is 



or 



1 + 77 TT^CUrlTT* ) d^lTri ( g ^ ) + 

J TT^CUrlTT*^ I j 



5 ° ( 7r "(o)!)<)^(^) + 

[I + 77 7r^curl7r* j 7r^curl7r* ( ^ / ^ ( 1/ ) = + ^ ® 



o o 

P^curl c curlP^ 



5o( q °) + (l + 77 curl) curl^ ^ ) ^ ^ = fl + 77 curl") (.7 + 5 ® W ) . 

Discussing this evolution equation we can finally make rigorous sense of our above argu- 
ments. 

The solution theory of the original Drude-Born-Fedorov model can now be formulated 
within the framework of the Sobolev lattice 
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2 Discussion of the Drude-Born-Fedorov Model 



H Stk [R,H S [ curl + i 

' k,. 

Theorem 2.1. Under Hypothesis Q and assuming that 



— e ( p ( curl ) n R ) u \ A| A isolated in a ( curl ) [• (12) 



'/ 

we have that for 

J e H e>0 \R,R[l + n curl ) ) , J = on R <0 , (13) 

and 

£/ie Drude-Born-Fedorov model 

^MclTK" 1 ^^ (15) 

' D\ A »,\ /e0\ 



b; V 1 + ' cur Vlo JU. 1 (16) 

/ias a unique solution (^^^j G Hq,o 0^> ^ 2 (^)) satisfying 

u 

in H-i ^curl + ij . TTie solution depends continuously and causal on the data. 
Proof. We have that 



' (0+) = Wo 



curl,, : D |^curlj r, R yl + 7] curlj c R M + 77 curlj —> i? ( 1 + 77 curl J c i? ( 1 + 77 curl 



h- > curl^ 

is a selfadjoint operator with a well-defined continuous inverse 



1 + 7] curl,, ) : R [ 1 + 77 curl ] — > ( 1 + 77 curl ] c i? ( 1 + 77 curl ] . 



It is 



7T^CUrl7T* = curl,, 



where ir, n denotes the canonical projection onto the subspace R I 1 + 77 curl ] of L 2 (Q) 



Then 7r* is the identity embedding of R yl + rj curlj into L 2 (O). 
We consider initially the solution of the evolution equation 



£ 



far, ( q " J + ^ 1 + r) curl,, J curl,, ^ J Q 



-1 



(17) 



1 + 77 curl,, J" + 



-i 



+5® y- + rj curl,, J W 

and observe that the resulting system is now an evolutionary system in the sense of Section 
[□with 



A = and 



C = (l + ?7 curl,, j curl,, ( ^ ^ 







I 1 + 77 curl,, I curl,, 



1 ° \ 
l + r/ curb, ) curb, 





continuous on R [ 1 + rj curl . 



I + 77 curl,, ) curb, = 77 I 1 + 77 curb, ) I r/curl^ + 1 J — 77 I 1 + 77 curl 



/ x- 1 / - 

— 7? _1 (I + 77 curl,, 



From the abstract theory we now find existence, uniqueness of a solution in the space 



H g fi ^R, i? ( 1 + 77 curl ) 
Moreover, 
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2 Discussion of the Drude-Born-Fedorov Model 



k v ([J ° J < (jj- x K>0 ® ( 1 + r/curl, ] W e H g>1 \ R,R[1 + rj curl 

With 

and using that P v = 7r,*7r^ commutes with ^ ^ \ we get 



Applying I 1 + curl now yields 



l + ,curl](^)(^] + 



-X R>0 (g> ^1 + V curl) 7T* ^1 + 77 curl^ W e H gA (r, #_i (curl + i) ) 

Clearly, 

1 + r] curl) 7r* ^l + r/ curl,) W = ix ri ^1 + rj curl) 7r* ^1 + rj curl,,) W 



1 + rj curb,) ^1 + r) curl,,) Wo 

= w 

and so 

l + ?7 curl) f £ Q °) (E\- XK>o ® Wo eH gA (r,H^ (cutrl + i 
The latter shows by causality and a temporal Sobolev embedding argument that 
t~((l + V cml) ( £ n °) (D-X^VWo) (t) 



ji J \H 



l + ?7 curl) ( e °) (l ) -v (x)F+ )(0) = 0. 



is continuous on IR and in particular 

fi J \H 
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The latter implies 



(0+) = W . 



Applying I 1 + 77 curl J 7r* to equation ( ITT)) yields similarly 



D 



-1 \ / E 



J + 5®W , 



which is the equation of the Drude-Born-Fedorov model. Note that this equation now holds 

o 

in if e __ 1 (R, iif_i(curl+i)). Causality and continuity estimates follow from the general theory 
applied to (TTTj) . 

Conversely, if ( ^ J solves ( TT5|) with ( TT6j) with J" = and Wo = 0, then 



do ( 1 + 77 curl 



e0\ * . / -1 

oJ +cml (i 



(19) 



and so 



curl) (^0 (q ° ) + ?7 ] 



Thus, we have that 



This observation implies 



1 



-i/O -1\ 



1 \H r 



1 7 V # 



e i? ( 1 + 77 curl . 



(20) 



1 



E 
H 



1 + 7] curl,, I T] 



-1 



1 U 



or 



( ^ ) 



-1 

1 



0. 



This implies 



= <He < R 



E 



( u ) + CUrl?? ( 1 + 77 CUllv 



-1 



1 



E 



n \ H 



» H 

E 
H 



ro 



e0\ E 
0» V \H 



5,0,0 



e0\ 

on) "U 



5,0,0 



(3,0,0 
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3 A Note on Generalizations of the Drude-Born-Fedorov Model 



and we see 



; h 1 



On the other hand, we read off from (1201) that 



= (1 - P TI ) 



-1\ ( E 

1 ) [ H 



implying 



E 
H 



This shows I 1 = 0, i.e. uniqueness. □ 



H 

Remark 2.2. Since the known results on the spectrum of curl show compact resolvent in the 



ortho-complement of the null space of curl, we also have R yi + r\ curlj = R yl + r/ curl 
in these cases. If we have — - e p I curl J the result becomes particularly simple. In this 



case we always have R ( 1 + rj curl^ = R \1 + r\ curl^ = L 2 (fi), which eliminates the 
possibly undesirable range conditions f lT3|) . fll4|) . 

Note that the assumption of Theorem 12. II may never hold (making the claim of the theorem 
trivially correct), as e.g. in the exterior domain case. However, for example in the case 
of Q having bounded measure, we have - according to [5] - pure point spectrum with no 
accumulation points and so the assumption of Theorem 12.11 is non-trivial and says simply 



-i e R or 

v 



e R\{0}. 



In the case 77 = 0, although not covered by the theorem, the solution theory is just the 
standard regular case. 



3 A Note on Generalizations of the Drude-Born-Fedorov 
Model 

By virtue of the power of the theoretical framework the above solution strategies can be 
conveniently generalized. We may straightforwardly generalize the Drude-Born-Fedorov 
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model ( I16p by imposing instead 

(b) = { k ^ 1 ) + c ^ 1 ) m *^(h) (21) 

where 
and 

with ^*,o selfadjoint, continuous, strictly positive definite and k, M* bounded, analytic 
in Be (r,r), for some r e lR>o- It should be noted that convolution integral type gener- 
alizations of the Drude-Born-Fedorov model have been previously suggested in |12) . but 
discussed only for the time-harmonic case. Also, the particular model based on homoge- 

nization with k (Sq 1 ) = r]^ 1 , M #i i = rje kx, M^o = ? ^^q^^' ^ e nas been suggested 



in [21]. 

Assuming 



as a replacement for assumption (lT2"j) we find that (<9 + curl J is boundedly invertible 
for all sufficiently large g e M > j^. Indeed, in this case 



■1 6 p ^k X curl^ 



with 



(Oq 1 ) + curl 
K (Oq 1 ) + curl 
= - <W (d Q l ) + curl^j 
= (l-Qo {d, 1 ))' 1 ^o + curl^ 

Qo (do 1 ) ■= (ko + curl^j d^Ki (Oq 1 ) 



as an ad-hoc abbreviation. Moreover, 

-i 



K (d x ) + curl J curl = 
(l - Qo (do 1 )) 1 (^-0 + curl^ curl 
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3 A Note on Generalizations of the Drude-Born-Fedorov Model 
shows that 



o \ 1 o 



C = in (d x ) + curl J curl 
is a bounded operator. So, with our assumption —lepl /tocurl ) , as a by-product of the 



proof of Theorem 12. 1\ the well-posedness result extends also to this generalization. The 
modified material law for the analogue of (|T7|) is 

M (d^ 1 ) = M* (d ') + d^C. 

Note that the initial condition now assumes the form 

(U ^o 1 ) + curl) M* (do 1 ) (Jj J J (0+) = L + curl) M 0) * (0+) = W . 

In summary, we obtain as a simple corollary to the proof of Theorem 12.11 (taking Remark 
12.21 into account): 

Corollary 3.1. Let Hypothesis fl hold and let 

k (Bq 1 ) = k + d^Ki (d^) 

and 

M* (d^) = M*, + do'M^ (d^ 1 ) 

with Kq, M^.o selfadjoint, continuous, strictly positive definite and k, M* bounded, analytic 
in Bq (r, r), for some r e M >0 . Then we have that for every 

- i e p ^curl) n R, (22) 



and 



i/ie Drude-Born-Fedorov model 



Jg^ (R,L 2 (Q)) (23) 
Wo e L 2 (fi) (24) 







*m+ ° " curl if ) = j+*®w„ (25) 



/ \ curl 

with 



5 



« (a - X ) + curl) M* (rV 1 ) ) (26) 
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References 



has a unique solution 



E 
H 



) 



eH e>0 (ffi, L 2 



(O)) satisfying 




o 



in H_i ( curl + i . Moreover, the solution depends continuously and causal on the right- 



hand side data. 
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